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General rogue waves in the Davey-Stewartson-II equation are derived by the bilinear method, and 
the solutions are given through determinants. It is shown that the simplest (fundamental) rogue 
waves are line rogue waves which arise from the constant background in a line profile and then 
retreat back to the constant background again. It is also shown that multi-rogue waves describe the 
interaction between several fundamental rogue waves, and higher-order rogue waves exhibit different 
dynamics (such as rising from the constant background but not retreating back to it). A remarkable 
feature of these rogue waves is that under certain parameter conditions, these rogue waves can 
blow up to infinity in finite time at isolated spatial points, i.e., exploding rogue waves exist in the 
Davey-Stewartson-II equation. 



^ I. INTRODUCTION 

o 

Rogue waves are large and spontaneous nonlinear waves and have been found in a variety of physical systems (such 
^ as the ocean and optical systems) [TJ |2] • Rogue waves generally occur due to modulation instability of monochromatic 
waves. One of the simplest mathematical models for modulation instability is the nonlinear Schrodinger (NLS) 
equation. For this equation, explicit expressions of rogue- wave solutions have been obtained by a variety of techniques 
^-H such as the Darboux transformation, the bilinear method and so on These NLS rogue waves can also be 

obtained from homoclinic solutions of the NLS equation under certain limits [T^HIS], or from rational solutions 

l_H of the Davey-Stewartson equation through dimension reductions [11] [17]. Physically these NLS rogue waves have 
C/^ been observed in optical fibers and water tanks [TH [19] . In addition to the NLS equation, rogue waves have also 
(— I been obtained in other wave equations, such as the Hirota equation, the derivative NLS equation and the Davey- 
• ^ Stewartson-I equation [5DH23]. Explicit rogue- wave solutions in mathematical model equations reveal the conditions 
^ for rogue- wave formation and facilitate the observation and prediction of rogue waves in physical systems [H [21 [TBI [E] ■ 
' In this article, we derive general rogue-wave solutions in the Davey-Stewartson-II equation. This equation arises 

^ in the modeling of two-dimensional shallow water waves Our derivation uses the bilinear method, and the 

^ solutions are expressed in terms of determinants. We show that the simplest (fundamental) rogue waves are line rogue 
waves which arise from the constant background in a line profile and then retreat back to the constant background 
l/^ again. We also show that the interaction between several fundamental rogue waves are described by multi-rogue- 
^-H wave solutions. However, higher-order rogue waves are found to exhibit different dynamics, such as rising from the 
constant background but not retreating back to it. An important feature about these rogue waves is that, under 
^Sj certain parameter conditions, these waves can blow up to infinity in finite time at isolated spatial points (we call such 
T-H solutions exploding rogue waves). The existence of exploding rogue waves is remarkable, and their appearance can 

be catastrophic in physical systems. 
' ^ [ It is noted that rogue waves are rational solutions of nonlinear systems in general. For the Davey-Stewartson 

^ equations, certain types of rational solutions have been derived before 1271 . Those rational solutions, under parameter 



X 



restrictions, would yield multi-rogue waves (see Sec. IIIB of this article). The rational solutions we would derive (in 



the next section), on the other hand, are more general; and these rational solutions, under parameter restrictions, 
could yield not only multi-rogue waves but also higher-order rogue waves. 



II. RATIONAL SOLUTIONS IN THE DAVEY-STEWARTSON-II EQUATION 



Evolution of a two-dimensional wavepacket on water of finite depth is governed by the Benney-Roskes-Davey- 
Stewartson equation |^iH25] . In the shallow- water (or long- wave) limit, this equation is integrable (see [15] and the 
references therein). This integrable equation is sometimes just called the Davey-Stewartson (DS) equation in the 
literature. The DS equation is divided into two types, DSI and DSII equations, depending on whether the surface 
tension is strong or weak |26j . 
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In this paper, we study the DSII equation. The normahzed form of this equation is 

iAt = A^.^ - Ayy + (e|A|2 - 2Q)A, 
Qxx Qyy — ^xx-i 
where e — 1 or —1. Through the variable transformation 

^ = ^y, Q = e-(21og/),„ 

where / is a real variable and g a complex one, this equation is transformed into the bilinear form, 

{Dl - Dl - iDt)g • / = 0, 
{Dl + Dl)f-f^2e[p-\g\^). 



(2.1) 



(2.2) 



(2.3) 



Rogue waves are rational solutions under certain parameter restrictions. Thus we first present general rational 
solutions to the DSII equation in the following theorem. The proof of this theorem is given in Appendix A. 



Theorem 1 The DSII equation (2.1) admits rational solutions (2.2) with / and g given by 2A^ x 2N determinants 

/ = To, 5 = (2.4) 

where 



(«) ^(n) 
i-n) 



(2.5) 



(«) 



fc=0 



(2.6) 



-An) 



fc=0 



/=0 



(2.7) 



Pi - £/Pi , Pt + e/Pi 
X H 



(2.8) 



^.7 



(2.9) 



the overbar ' ' represents complex conjugation, i,j = 1,...,A^, ni^mj are arbitrary non-negative integers, and 
Pi, qj, Cik, dji are arbitrary complex constants. 

Remark 1. By a scaling of / and g, we can normalize Cio = djo — 1 without loss of generality, thus hereafter we 
set do = djQ = 1. 



Remark 2. For e = —1, / in (2.4) is non-negative, i.e., / > 0. A proof is given in Appendix B. Since / is the 



denominator of the solutions A and Q, the above rational solutions are nonsingular as long as / > 0. But it is also 
possible that / hits zero and the corresponding solution blows up to infinity at a certain point of space-time, which 
we will see later. 

Remark 3. Rational solutions in the DS equations have been derived in [37] before. The nonsingular rational 
solutions for the DSII equation in that paper correspond to special rational solutions in the above theorem with 
ni = • • • = riN = 1 and mi = • • • = TOjv = 0. 

The simplest rational solution is obtained when = 1, ni = 1 and mi = 0. In this case. 
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where 



An) 



Pi + qi 



pi 



Pi +qi 



Cu 



An) 



1 , piqi 

,4i +"-- ^ , _ +Cii), 



+ e 



and is defined in (2.8). This solution seems to have three free complex parameters pi, qi and cn, but qi can be 
absorbed into cn by a reparametrization. Indeed, by defining 



= cii 



Pi 



(biP"e)(kip-e) V Pi +11 



Piqi + e 



and denoting pi = p, Ci + = S,, we can show that the terms in t„ which are linear in ^ + n and £, — n vanish, and 
this Tn reduces to 

Tn = (e + A, 



^ = ax + by + ujt + e, A 



-e\p\' 



p—e/p , p + e/p. p^ + l/p^ 
a=^^, h=^^^, uj= , , 

up to a constant multiplication, thus this new r„ yields the same solution. The solution from this new t„ has only 
two independent complex parameters p and B now. If we separate the real and imaginary parts of a,b^uj and 9 as 

a = fli + ia2, h = bi + ib2, uj = Ui + iuj2, 9 ^ 9i + i92, 

then the explicit expressions for this solution are 

2i{a2X + b2y + i02t + 02) + 1 



A{x,y,t) = V2 



1 



/ 



(2.10) 



where 



Q(x,y,i) = e-(21og/), 



/ = {aix + biy + Wit + 9if + {a2X + b2y + + ^a)^ + A. 



(2.11) 



This simplest rational solution is nonsingular when e = — 1 (where A > 0). In this case, the solution exhibits two 
distinctly different dynamics depending on the parameter value of p. 

1. If \p\ 7^ 1, then it is easy to see that b/a is not real, hence 61/62 7^ 01/02. In this case, along the [x{t),y{t)] 
trajectory where 

oix + 6iy = -wii, a2X + b2y ^ ~uj2t, 

solutions {A,Q) are constants. In addition, at any given time, {A,Q) — {^/2,e) when {x,y) goes to infinity. 
Thus the solution is a two-dimensional lump moving on a constant background [27 . 

2. If \p\ — 1, then a,b are real but w is imaginary. In this case, the solution depends on {x,y) through the 
combination aia; + biy and is thus a line wave. As t — ^ ±00, this line wave goes to a uniform constant 
background (as long as p^ ^ ±z); in the intermediate times, it rises to a higher amplitude. Thus this line wave 
is a line rogue wave which "appears from nowhere and disappears with no trace" . 



When e = 1, the rational solution ( 2.10 )-( 2.11) is singular on a certain elliptic curve in the {x,y) plane for any time 
t, since A < now. For this e, the constant-background solution is modulationally stable [22] i thus no rogue waves 
can be expected. In view of this, we only consider the case of e = — 1 in the remainder of the paper. 
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III. ROGUE WAVES IN THE DAVEY-STEWARTSON-II EQUATION 

As we see from the above analysis, rogue waves would result from the rational solutions in Theorem 1 for e = — 1 
under certain parameter conditions. Specifically, to obtain rogue waves, we need to require e = — 1 and 



1, if Uj > 0; 



1, if mj > 0; 1 < j < A^. 



(3.1) 



In this section, we examine the dynamics of these rogue waves in detail. 



A. Fundamental rogue waves 



Fundamental rogue waves in the DSII equation are obtained when one takes 

e = -1, = 1, ni = 1, mi = 0, pi = e''^ 



(3.2) 



in the rational solution (2.4), with j3 being a real parameter and pf ^ zti (i.e., cos2/3 7^ 0). As we have explained in 
the previous section, this solution is equivalent to ( 2.10| -(2.11 1. After a shift of time and space coordinates, 9i and 
02 can be eliminated. Then in view oi p — e*^, this fundamental rogue wave becomes 



A{x,y,t) = V2 ( 1 



4- 16itcos2/3 



1 + 4(x cos /3 - y sin /?)2 + IGt^ cos^ 2^ 



(3.3) 



Q{x,y,t) = -1 - 16cos^/3 



1 - 4(x cos /? - 2/ sin /3)2 + 16^^ cos^ 2/3 
[1 + A{x cos /3 - ?/ sin /3)2 + 16^2 cos2 2/3]^ 



(3.4) 



where /3 is a free real parameter. This solution describes a line wave with the line oriented in the (sin /3, cos /3) 
direction of the {x,y) plane, and the orientation angle is 7r/2 — /3. The width of this line wave is the same for 
all /? values, i.e., the width is angle-independent. At any given time, this solution is a constant along the line 
direction (with fixed x cos P — y sin /?) and approaches the constant background away from the center of the line (with 
xcosP — ysin/3 — )■ ±cx)). When t — )■ ±c», the solution A uniformly approaches the constant background v^; but 
in the intermediate times, |A| reaches maximum amplitude 2>^/2 (i.e., three times the background amplitude) at the 
center [xcosfi — ysin/3 = 0) of the line wave at time t = 0. The speed at which this line wave climbs to its peak 
amplitude is proportional to | cos2/?|, which is angle-dependent. This fundamental rogue wave is illustrated in Fig. 1 
with /3 = 7r/6. 

It is noted that under the same parameter conditions (3.2) but with cos 2/3 = 0, i.e., this line wave is oriented 
diagonally (45°) or anti-diagonally (—45°), then a; = in the rational solution (2.10)-(2.11). In this case, after a shift 
of space coordinates, 9i can be eliminated. Hence this rational solution becomes 

/ /I I \ 

(3.5) 



A{x,y,t)^V2[l- - 



2{x ± yf + 46*2 



Q{x,y,t) 



l~2{x±yf 



[l + 2{x±yy 



Aeir 



(3.6) 



where 62 is a free real parameter. This solution is not a rogue wave. Instead, it is a stationary line soliton sitting on 
the constant background. Its peak \A\ amplitude is y^2]^~+4^|)7(T-i-4^|y. The highest value of this peak amplitude 
is 3-\/2 (three times the constant background), which is attained at 62 = 0. When 16*21 increases to infinity, this peak 
amplitude decreases to the background amplitude ^/2. 

li N > \, ov N ^ 1 but mi -I- ni > 1, the rational solutions in Theorem 1 under parameter restriction ( 3.1 1 will give 
a wide variety of non-fundamental rogue waves. For simplicity, we consider three subclasses of such solutions below. 



B. Multi- rogue waves 

One subclass of non-fundamental rogue waves is the multi-rogue waves which describe the interaction between 
several fundamental rogue waves. These solutions can be obtained from Theorem 1 by taking 



-1, 7V>1, 



1, 



0, p, 



1<J<N, 



(3.7) 



< 2 



< 2- 



i = -10 




t = 




< 2- 



-10 



< 2 



^ = -1 



10 



i = 10 




FIG. 1: A fundamental rogue wave (|3.3[) with /3 = 7r/6 



where /3j is a free real parameter (with cos2/?j 7^ 0). In this case, the r-solution (2.5) becomes 



where 



J") 



1 



m. 



(«) 



C + " - 



Pi + (13 



(3.8) 



(3.9) 



(3.10) 



is defined in (2.8), and Qj, ca are free complex constants (but with qj 7^ ±pi to avoid zero divisors). When t — )■ ±cx), 
the solutions (A7Q) approach the constant background uniformly in the entire (x, y) plane. In the intermediate times, 
N fundamental line rogue waves arise from the constant background, interact with each other, and then disappear 
into the background again. Depending on the parameter choices, individual line rogue waves can reach their peak 
amplitudes at the same time or at different times, with the former yielding stronger interactions. 

Now we illustrate these multi-rogue waves and examine their dynamics. To obtain a two-rogue wave solution, we 
take parameter values 



= 2, pi = 1, _P2 = «, <?! =0, 52 = -3, cu = 0, C21 = i/2. 



(3.11) 



The corresponding solution \A\ is displayed in Fig. [2] It is seen that as i ^ ±cx), the solution uniformly approaches 
the constant background v^; but in the intermediate times, a cross-shape rogue wave appears. This cross rogue wave 
describes the interaction between two fundamental line rogue waves, one oriented along the y direction (corresponding 
to the parameter pi), and the other one oriented along the x direction (corresponding to the parameter y^)- These two 
individual line waves reach their peak amplitude 3-\/2 at different times, with the x-direction one peaking at t w — 1/4 
and the ^/-direction one peaking at t w 0. 
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t=-4: 



t = -0.8 



t= 0.35 




10 -10 



10 -10 



10 -10 



FIG. 2: A two-rogue wave solution (3.81 with parameters (3.111 



Next, we take parameter values 

iV = 4, pi = 1, p2 



-0.1, 92 = 0, 



(3.12) 



93 = 0.1, 94 = 0.2, cii = -2i, C21 = 0, C31 = 2i, C41 = 



(3.13) 



which gives a four-rogue wave solution. This solution (|^|) is displayed in Fig. 3. As t — > ±cxi, the solution uniformly 
goes to the constant background ^/2] but in the intermediate times, a rogue wave comprising four lines emerges. 
These four individual line waves reach their peak amplitudes 3-^/2 at approximately the same time i = 0, and their 
widths are identical (see t = Q panel). Due to the interaction of these four line waves, the maximum amplitude of the 
solution (at intersections of the four lines) can be very high. Indeed, at t = — 1, we find that the peak amplitude of 
the solution |^| reaches approximately 30\/2 (i.e., 30 times the constant background). Thus such rogue waves can be 
fairly dangerous if they arise in physical situations. 

In the general iV-rogue wave solution (3.8), /3j is a free real parameter, and qj,Cji are free complex parameters 



(1 < j < N). Thus it appears that this A^-rogue-wave solution contains A^ free real parameters and 2N free complex 
parameters, totaling 5A^ free real parameters. But these parameters are reducible (similar to the simplest rational 
solutions in the previous section). Indeed, when A^ = 2, by a reparametrization of 



Oil = Cii - p, 



{i:( 



1 



1i 



1 



1 



Pi - P3- 



i = l,2. 



we can show that the two-rogue-wave solution (3.8) is reduced to 



((0 + + n) - ((Ci - ..)(& - ») - j;^) 



1 1 

2bi+P2| 



1 

16 



Pi -P2 
Pi +P2 



(3.14) 



up to a constant multiplication (which does not affect the solution). Here d = £,[ -\-Cii, z = 1,2. In this equivalent 
T„ solution, parameters qj disappear, thus it contains only /3i, ^2, cn and C21. Of the two complex constants cn 
and C21, their real parts and the imaginary part of one of them can be further normalized to be zero by a shift of 
the {x,y,t) axes. Thus this two-rogue wave solution contains only three irreducible real parameters. For the general 
A'^-rogue wave solution (3.8), we conjecture that all qj parameters can also be removed by a reparametrization of Cn, 
hence this A^-rogue wave solution contains only 3(A'^ — 1) irreducible real parameters (after a shift of {x,y,t)). 



y 




FIG. 3: A four-rogue wave solution ( 3.8 1 with parameters ( 3.12 l-( 3.13 1. Plotted is the \ A\ field, with red color indicating higher 
values. The spatial region in each panel is — 20 < a::, y < 20, and the constant-background value is ^/2. 



It is noted that if instead of (3.7), one takes 



> 1, rii = 71,2 



TT-AT = 0, mi = m2 



ruN = 1, 



then the same multi-rogue-wave solutions as above will be obtained. Thus different parameter choices can lead to the 
same solutions. 

C. Higher-order rogue waves 

A second subclass of non-fundamental rogue waves is the higher-order rogue waves. These solutions are obtained 
from Theorem 1 by taking 



-1, iV = l, ni > 1, TOi=0, bi| = l. 



(3.15) 



In this case, the r-solution (2.51 becomes 



(«) 

™ii rrii^ 



^(-n) (-«) 



(3.16) 



where 



"1 H 



fc=0 



Pi + qi 



(3.17) 



= > cifc(pi5p, +^[ +n) 



ni—k 



1 



k=0 



(3.18) 



(,1 is defined in (2.8), cio = 1, and cik,qi are free complex constants. These higher-order rogue waves exhibit dynamics 



different from those of multi-rogue waves, as we will demonstrate below. 



8 



For simplicity, we consider second-order rogue waves where ni = 2. In this case, we find that 



„2 



2 



m 



(n) 



Pi9i 



piqi - 1 



Cll 



+ Cl2 



-11 



+ 



Pi9i 



{PiQi - 1) 



1)2}' 



where 



Denoting 



Ci ^Pi^pifi = ^ X 



Pi + iM ./— y I Pi- i/Pi 



2t. 



P = Pi, 



91, e = + C = ei 



where a = Cii/2 — 1 and 6 = C12 — the r„ sohition (3.16) becomes 



1 



p + q 
1 n2 



pq-l 



+ C + 



(p + g)2 
pq 



}{(^ 



9 



P9 



1 



+ C + 



{p + g)2 
pq 



{pq - 1) 



1)2}- 



This solution has four apparent complex parameters, p, g, a and b. But q can be removed by a reparametrization of a 
and 6. Indeed, by replacing 



a—l+p 



1 



b~¥b + p 



g P 

— 1 |p|2 + 1 
9 P 



,(p + g)2 (to -1)2 (bP + l)2. 
and recalling \p\ — 1, the above r„ can be rewritten as 

je + ")' + c)((e-n)' + c) + + (3.19) 

up to a constant multiplication. Thus this second-order solution contains only parameters p, a and b now. 

In these second-order solutions, if p^ ^ ±i, then the solutions do not uniformly approach the constant background 
as < — ±00, thus they are not rogue waves. But when p^ = —i, the solution uniformly approaches the constant 
background as i — >■ —00, thus it "appears from nowhere" and is a rogue wave. However, this second-order rogue 
wave does not retreat back to the constant background when t — >■ -l-oo, thus it does not "disappear with no trace". 
This means that this second-order rogue wave behaves quite differently from the multi-rogue waves considered in the 
previous subsection. 

Below we examine this second-order rogue wave in more detail. For definiteness, we take p = e^*'^/^ (the choice of 
p — —e"™/'^ would yield the same solution). In this case, by a shift of {x, y, t) axes, we can normalize b as well as the 
real part of a to be zero. Thus we can set 



a = ia, b = 0, 



(3.20) 



where a is a free real parameter. Substituting these p, a and b values into the t„ solution (3.191, we find that the 
solution A{x,y,t) becomes 



A = V2 



(1 + 2ia)[{x + yf + 8t] - 2i{x'^ - y'^ + a - 2a^) + 6a'^ 



At - a 



+2(a(2; + y) - i(a;-y)) + l{x + y)^ + 



(3.21) 



and the solution Q{x, y, t) is given by (2.2 1 with e — —1 and / being the denominator in the above A solution. When 
t — > —00, this solution A{x,y,t) uniformly approaches the constant background \/2 (like regular rogue waves). But 
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when t — )■ +c», it approaches two lumps which slowly move away from each other. The peak amplitudes of these two 



lumps are attained at {x,y) locations where the first two terms in the denominator of (3.211 vanish, i.e., at 



a^max 



= ±(i +a) Vst + 2a^, 



J/max 



and these peak \A\ amplitudes approach 3-\/2 when t — > +00. This solution with a = 1 is displayed in Fig. [4j We 
see that this second-order rogue wave looks quite different from the previous rogue waves in Figs. 1-3. Instead of 
"disappearing with no trace" , this second-order rogue wave "disappears with a trace" . 



t= 20 



t=-l 



t = 



t = 2 



t = -0.b 




t = 10 




FIG. 4: A second-order rogue wave solution (3.21 1 with a = 1 



It is noted that when — i, i.e., p = zte"^/"*, the second-order solution (3.19) would approach the constant 

-c». In other words. 



background when t — > -l-oo but approach two lumps which move away from each other when t 
this solution describes a process which is opposite of that when — —i (see Fig. m. 



D. Exploding rogue waves 

A third but important subclass of non-fundamental rogue waves is the exploding rogue waves. These rogue waves, 
which arise from the constant background, can blow up to infinity in finite time at isolated spatial locations. These ex- 
ploding rogue waves can be obtained from the higher-order rogue waves or multi-rogue waves under certain parameter 
conditions, as we will demonstrate below. 



First, we consider the second-order rogue waves (3.21). When a = 0, this solution becomes 



A{x,y,t) = V2 



1 - 



{x + y)'^ + 8t - 2i{x^ - y2) 



(3.22) 



This solution uniformly approaches the constant background \/2 as i — > ±00. But in the intermediate time t — 0, it 
blows up to infinity at the origin {x, y) — (0, 0). To see this, we notice that at (x, y) = (0, 0), 



A(0,0,i) = 



2t 



(3.23) 



thus this A solution blows up to infinity when t approaches zero (the solution Q blows up to infinity at this time as 
well). The rate of blowup is (t — t^,)~^, where = is the time of singularity. This exploding process is displayed in 
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i = 10 



i = 0.5 





t= Q.l 



i= 0.01 




FIG. 5: An exploding second-order rogue wave (3.221 



Fig. [5] The existence of exploding rogue waves in the DSII equation is a remarkable phenomenon, and their occurrence 
would be catastrophic in physical systems. 

In addition to higher-order rogue waves, multi-rogue waves can also explode under suitable choices of parameters. 
To demonstrate, we consider the two-rogue- wave solutions whose simplified expressions are given in Eq. (3.14). Taking 
parameter values 



Pi = 1, P2 



cii = C21 = 0, 



(3.24) 



this two-rogue wave becomes 



where 



A{x,y,t)^V2-, Q = -l-(21ogTo), 

To 



(3.25) 



2 2 
Ti= X y 



At the origin (x, y) — (0, 0), 



+ + 4^2 + ^ ) - Ait{x^ - 2/2). 



A(0,0,t) = \/2 



i2 



(3.26) 



thus this wave explodes to infinity at times = ±-\/3/4. Its exploding rate is also (t — t*)^^, where is the time of 
wave singularity. This exploding two-rogue-wave solution is displayed in Fig. [6j 

It is noted that for the Davey-Stewartson equations, self-similar collapsing solutions have been derived in [501151] . For 
the non-integrable Benney-Roskes-Davey-Stewartson equations, wave collapse has also been reported 1351 [33]. Those 
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i = 10 



t= 1 





t= Q.6 



t = 0.44 




FIG. 6: An exploding two-rogue wave (3.251 



collapsing solutions are different from our exploding rogue waves since the boundary conditions of those solutions are 
different. Specifically, those collapsing solutions do not arise from the constant background and are thus not rogue 
waves. 

In this section, only a few subclasses of rogue- wave solutions were examined. The rational solutions in Theorem 1, 
under parameter conditions (3.11, also contain a lot of other subclasses of rogue waves which are not elaborated in 



this article. We also note that different choices of parameters can yield the same solutions. For instance, if we take 

A^ = l, ni = mi = 1, \pi\ = \qi\ = l 



in Theorem 1, the resulting solution (2.4) would be equivalent to the two-rogue-wave solution (3.8) with parameters 

N = 2, 712 = 1, mi = TO2 = 0, |pi| = |p2| = 1 



(see also Eq. (3.14)) 



IV. SUMMARY AND DISCUSSIONS 



In this article, we have derived general rogue waves in the Davey-Stewartson-II equation. We have shown that the 
fundamental rogue waves are line rogue waves which arise from the constant background in a line profile and then 
retreat back to the constant background again. We have also shown that multi-rogue waves describe the interaction 
between several fundamental rogue waves, and higher-order rogue waves exhibit different dynamics (such as rising 
from the constant background but not retreating back to it). In addition, we have discovered exploding rogue waves, 
which arise from the constant background but blow up to infinity in finite time at isolated spatial points. 

It is helpful to compare these rogue waves in the Davey-Stewartson-II equation with those in the Davey-Stewartson-I 
equation (see [13]). The biggest difference is that exploding rogue waves exist in the Davey-Stewartson-II equation, 
but such waves cannot be found in the Davey-Stewartson-I equation |23| . In Appendix C, nonsingularity of rogue 
waves in the Davey-Stewartson-I equation is analytically proved for a subclass of parameter values, and we conjecture 
that all rogue waves (which arise from the constant background) are nonsingular in the Davey-Stewartson-I equation. 
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Other differences on rogue waves also exist between the Davey-Stewartson-I and Davey-Stewartson-II equations. 
For instance, in the Davey-Stewartson-I equation, fundamental (line) rogue waves can only be oriented along a half of 
all possible angles in the {x, y) plane [53]; but in the Davey-Stewartson-II equation, fundamental rogue waves can be 
oriented along any angle (except diagonal and anti-diagonal angles). This difference has important implications for 
multi-rogue-wave patterns. For instance, in the Davey-Stewartson-II equation, cross rogue-wave patterns formed by 
two orthogonally-oriented fundamental rogue waves exist (see Fig. [2|; but in the Davey-Stewartson-I equation, cross 
patterns of multi-rogue waves cannot exist. 
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Appendix A 



In this appendix, we derive the rational solutions to the DSII equation in Theorem 1. 
The bilinear form (2.3) of the DSII equation can be derived from 



{Dx^Dx_i — 2)t„ • T„ = — 2t„+it„_i, 

{Dl^ - D^^)Tn+l ■ Tn = 0, 



(A.l) 



2 

X-1 



Dx_2)'Tn+l • T„ = 0, 



by taking the independent and dependent variables as 



Xl 



2i 



t, X_2 = 



2i 



(A.2) 



(A.3) 



and imposing the complex conjugate condition 



The variable transformation (A.2) means that 



x_i = exi, X_2 = X2. 



(A.4) 



(A.5) 



We first consider the ratio nal solutions for the system (A.l ), and then obtain those for (2.3) by imposing the complex 
conjugate condition (A.4). 



It is known that the bilinear equation (A.l) admits determinant solutions 



det 



(A.6) 



where K is a positive integer, toIJ"* is an arbitrary function satisfying the differential and difference relations. 



(A.7) 
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and <^|"\ V'l"'' arbitrary functions satisfying 



o (n) (n+i') o ; (n) , (n— i^) / i i i ri\ 

^V'j , - -V'} , (J^ = ±l,±2). 



The rational solutions are obtained by taking 



1 1 2 

Pi Pi 



(A.8) 

(A.9) 

(A.IO) 
(A.U) 



1 1 2 

??j = -^a;_2 H x_i + qjxi - q.X2, 



(A.12) 



where pi and qj are complex constants, and Bj are differential operators of order rii and with respect to pi and 
qj respectively, defined as 

A,; = X] - ^ (A.13) 



fe=0 



Cjfe, c?j! are complex constants, and n^, are non-negative integers. It is easy to see that the above 'mf^\ tpf"^ and 



ip'f^^ satisfy the differential and difference relations (|A.7|) and (|A.8| 



Next we impose the complex conjugate condition (A. 4) wit 



general rational solutions (A.6I with 2N x 2N determinants (i.e., K = 2N) 

T„ = det ( to|"^ ) , 



1 the restriction (A.5). For this purpose, we consider 

(A.14) 



together with (A.9) and (A.13). In this solution, we impose the parameter conditions 



e e 

PN+i — — , iN+j = — , n-M+i — rii, rriN+j — mj, 

Pi qj 



k 

cN+i,k = 



for 1 < i-, j < N . Under these parameter conditions, we have 

^N+i = ^i, rjN+j = fjj, 



fe=0 Ai=0 



1=0 



u=0 



Using the operator identities 



fe- dp, -I) ^Pj = Pj ipj 9p, ) ^ {qj 9g, - 1 ) = 9j ) , 
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the elements of the determinant in r„ become 



An) 



N+i ~ ^i^N+j 



9j 



(") _ 

(«) 
''■N+i,N+j 



1 



1 



Pi 



i+Vj 



Pi + V qj . 

Since the t„ solution can be scaled by an arbitrary constant, we define a scaled r„ function as 

N 



rn/l[{-ir+"''ep,q,^T^. 

i=l 



This scaled r„ solution can be written as 

(n) 

T, 



(-1)"^' (") 



(-1)"' („) (-1)"'+"^-- 



i,N+j 



in) 
'^N+i,N+j 



(«) ^(n) 



where 



(-ir^-_(„) 



_(_„) (_„) 



(A.15) 



:{-mqjTe 



(A.16) 



We can s ee fro m (A.15 1 that this t„ satisfies the complex conjugate condition (A. 4 1, and thus it satisfies the bilinear 



equation (2.3) of the DSII equation. 

Finally we simplify the above r„ solution. Using the operator identities 

{PidpMe^^ = Pre«-(p,9p, + e^ + n), 



where 



fe9,,)(-g,)-"e"^- = (^q^)--e^'{q,d,^+v',-nl 
2 1 2 1 

^x_2 x^i +PiXi +2pfx2, Vj = ^a;_2 a^-i + QjXi - 2q'^X2, 

Pi Pi 1j 1j 



the rational solutions to the DSII equation can be obtained from (A. 13 1, (A.15 1 and (A. 16 1 as 



Tn = 



(n) ^(n) 



(A.17) 



where 



m 



in) 



7Li •"'J 

9i , ^ 



fc=0 



1=0 



Pi + 



An) 



77, + n 



1 



fe=0 (=0 

Then using the gauge invariance of r„ , we see that Tn with matrix elements 



Piqj + e 



niA—l 



k=0 



1=0 



Pi + Qj 
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m. 



(n) 



fc=0 



also satisfies the bilinear equation ( A.l ) as well as the complex conjugate condition (A.4|, thus it satisfies the bilinear 



equation (2.3) of the DSII equation. This completes the proof of Theorem 1. 



Appendix B 



In this appendix, we prove that / in Theorem 1 is non-negative for e = —1. In view of Eqs. (2.4| and (2.5), it 
suffices to show the following lemma. 



Lemma 1. For any N x N matrices A and B, the following 2A^ x 2A^ determinant is non-negative: 

> 0. 



A_ B 
-B A 



Proof. Wc will prove this lemma by induction. For = 1 the statement in the lemma is obviously true. Let us 
denote N x M matrices as 



i-NM 



mat 

l<j<Ar, l<j<M 



(ay) , 



B 



NM 



mat 

l<i<N. l<j<M 



where aij and 6y are complex numbers. Then by the Jacobi formula for determinants, we have 



B 



N+1,N+1 



— Sw+l.Af+l An+i^n+1 

An+i.n+i Bn+i^n 
Ann 

An+i.n Bn+i^n+i 



A 



NN 



B 



NN 



—Bnn Ann 
Ann Bn,n+i 

— Bn+1.N ^Af+l,7V+l 

An.n+1 Bnn 



-B 



NN 



A 



N,N+1 



-B 



N+l.N+l An+1^N 



(A.18) 



The right-hand side of this equation can be rewritten as 

2 



A 



N+1,N+1 Bn+1,N 



-B 



N,N+1 



A 



NN 



Ar 



N+1,N Bn+1,N+1 



-B 



NN 



A 



N,N+1 



which is non-negative. Denoting 



Dn = 



A 



NN 



B 



NN 



-B 



NN 



A 



NN 



then Eq. (A.18) gives Dn+iDn > 0. Therefore if Dn > 0, we get Dn+i > 0. If Dn — 0, then by an infinitesimal 
deformation (for example, Ann ~^ Ann + o.In with an infinitesimal real number a and the N x N unit matrix In), 
the deformed Dn becomes positive. Thus the infinitesimally deformed Dn+i is non-negative, and so is Dn+i- This 
completes the induction and Lemma 1 is proved. 



Appendix C 



In this appendix, we comment on the nonsingularity of rational solutions for the DSI equation given in [53^. The 
solution in Theorem 1 of is nonsingular if the real parts of wave numbers (1 < i < N) are all positive. Because 
if Re Pi > 0, then from the appendix in [23], it is easy to see that the denominator / is given by the determinant of a 
Hermite matrix whose element can be written as an integral, 



/ = det 

l<ij<Ar 



A^AjC 
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Here the condition of Repi > (for all 1 < z < A'') is used to guarantee that the antiderivative of e^'+^J (with respect 
to Xi) vanishes at xi = —oo. Then for any non-zero vector v = (ui, V2, - ■ ■ , vn) and *v being its complex transpose, 
we have 



^ 2 



i=l 



dxi > 0. 



This shows that the Hermite matrix ^m^^'^ is positive definite, hence its determinant / is positive, i.e., / > 0. 

When the real parts of wave numbers pi (1 < « < N) are all negative, by slightly modifying the above argument, 
we can show that the rational solutions in the DSI equation are nonsingular as well. 

We conjecture that the rational solutions in the DSI equation, as given in [23] , are actually nonsingular for all wave 
numbers pi (1 < i < N). 
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